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The NASA Short-term Prediction Research and Transition (SPoRT) Center is producing real-time, 1- 
km resolution Normalized Difference Vegetation Index (NDVI) gridded composites over a Continental 
U.S. domain. These composites are updated daily based on swath data from the Moderate Resolution 
Imaging Spectroradiometer (MODIS) sensor aboard the polar orbiting NASA Aqua and Terra satellites, 
with a product time lag of about one day. A simple time-weighting algorithm is applied to the NDVI 
swath data that queries the previous 20 days of data to ensure a continuous grid of data populated at all 
pixels. The daily composites exhibited good continuity both spatially and temporally during June and 
July 2010. The composites also nicely depicted high greenness anomalies that resulted from significant 
rainfall over southwestern Texas, Mexico, and New Mexico during July due to early-season tropical 
cyclone activity. 

The SPoRT Center is in the process of computing greenness vegetation fraction (GVF) composites 
from the MODIS NDVI data at the same spatial and temporal resolution for use in the NASA Land 
Information System (LIS). The new daily GVF dataset would replace the monthly climatological GVF 
database (based on Advanced Very High Resolution Radiometer [AVHRR] observations from 1992-93) 
currently available to the Noah land surface model (LSM) in both LIS and the public version of the 
Weather Research and Forecasting (WRF) model. The much higher spatial resolution (1 km versus 0.15 
degree) and daily updates based on real-time satellite observations have the capability to greatly 
improve the simulation of the surface energy budget in the Noah LSM within LIS and WRF. Once code is 
developed in LIS to incorporate the daily updated GVFs, the SPoRT Center will conduct simulation 
sensitivity experiments to quantify the impacts and improvements realized by the MODIS real-time GVF 
data. This presentation will describe the methodology used to develop the 1-km MODIS NDVI 
composites and show sample output from summer 2010, compare the MODIS GVF data to the AVHRR 
monthly climatology, and illustrate the sensitivity of the Noah LSM within LIS and/or the coupled 
LIS/WRF system to the new MODIS GVF dataset. 
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1. INTRODUCTION 

An important variable that especially impacts the 
transport of moisture into the atmosphere during the 
warm season is evapotranspiration from vegetative 
surfaces. Vegetation is represented in models by the 
horizontal and vertical distribution of plant vegetation 
given by the Greenness Vegetation Fraction (GVF) and 
Leaf Area Index (LAI), respectively (Gutman and 
Ignatov 1998). The operational Noah land surface 
model (LSM; Chen and Dudhia 2001; Ek et al. 2003) as 
found in the Weather Research and Forecasting (WRF) 
model (Skamarock et al. 2008) and National Centers 
for Environmental Prediction (NCEP) North American 
Mesoscale model (Janjic et al. 2001; Janjic 2003) holds 
the LAI fixed for all vegetation classes. 

The GVF, meanwhile, is allowed to vary spatially in 
the Noah LSM according to a global monthly 
climatology dataset derived from Normalized 
Difference Vegetation Index (NDVI) data on the NOAA 
Advanced Very High Resolution Radiometer (AVHRR) 
polar orbiting satellite, using information from 1985 to 
1991 (Gutman and Ignatov 1998; Jiang et al. 2010). 
Representing data at the mid-point of every month 
(e.g. June climatology is valid for the 15th of the 
month), this monthly climatological dataset is on a grid 
with 0.125° (~14 km) spatial resolution and has been 
implemented into the operational Noah LSM at NCEP 
and within the community WRF model (Ek et al. 2003; 
Jiang et al. 2010; Skamarock et al. 2008). 

A limitation that the climatological dataset 
presents is that the annual cycle of GVF is always 
represented the same in models from one year to the 
next. In reality, the response of vegetation to 
meteorological and climate conditions varies between 
seasons and years based on anomalous weather and 
climate features. Extreme events such as an unusual 
hard freeze, late bloom due to colder than average 
temperatures, or drought can lead to a vegetative 
response that is quite different than the climatological 
representation. In addition, the dated nature of the 
GVF climatology and relatively coarse resolution may 
not be representative of current vegetative conditions 
in today's high-resolution numerical models. Recent 
land use changes due to urbanization since the period 
of record of the GVF climatology likely contribute to 
mis-representations in the models. 
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Therefore, the NASA Short-term Prediction 
Research and Transition (SPoRT) Center proposes a 
Continental U.S. (CONUS) scale, high resolution GVF 
dataset that is updated on a daily basis with near real- 
time swath data from the Moderate Resolution 
Imaging Spectroradiometer (MODIS) instruments 
aboard the NASA Earth Observing System Aqua and 
Terra satellites. This new dataset is inserted into the 
NASA Land Information System (LIS) to improve the 
representation of land surface processes within the 
Noah LSM in LIS, and ultimately in the WRF model via 
the coupling to LIS. The GVFs are derived from NDVI 
data that are produced in near real-time from the 
Aqua and Terra platforms. The NDVI is based on 
properties of healthy vegetation, which has a high 
absorbance (low reflectance) in the visible portion (or 
photosynthetically active region) of the 
electromagnetic spectrum while having a high 
reflectance at the near-IR wavelengths. Thus, the NDVI 
is defined as the combination of these reflectances: 

NDVI = Pnir ~ Pred (1) ? 

Pnir + Pred 

where p N , R is the satellite reflectance at near-IR 
wavelengths (0.75-1.5 jum) while p RED is the satellite 
reflectance at visible-red wavelengths (0. 6-0.7 jum). 
NDVI ranges from -1 to +1, where a value near +1 
indicates a healthy, fully-vegetated surface. Values 
near 0 indicate little to no vegetation while negative 
values typically correspond to snow or ice cover under 
clear sky conditions. 

Experiments with real-data and/or near real-time 
vegetation in place of the operational climatological 
dataset are not new. Previous studies have examined 
near real-time vegetation datasets derived from the 
NOAA/AVHRR satellite (Jiang et al. 2010) and its 
potential utility on real-time modeling (Crawford et al. 
2001; Kurkowski et al. 2003; James et al. 2009). Few 
studies, however, have examined such datasets 
derived from the NASA MODIS instruments (Miller et 
al. 2006; Ruhge and Barlage 2011), particularly with 
real-time MODIS data at high spatial resolution (1-km) 
for regional, real-time high-resolution modeling 
applications (e.g. Kain et al. 2010). This paper and 
companion presentation describes the development of 
a CONUS-scale, 1-km real-time MODIS NDVI/GVF 
product to replace the coarser-resolution 
climatological GVFs currently implemented in the 
NASA LIS and operational models. Section 2 gives 
some background information on the NASA SPoRT 
Center. Section 3 provides a brief overview of the 
NASA LIS land surface modeling framework. Section 4 



describes the real-time production of the SPoRT NDVI 
composites and how the GVF is computed for use 
within LIS. Preliminary results from offline LIS-Noah 
runs are presented in Section 5 followed by a summary 
and future direction presented in Section 6. 

2. NASA SPORT CENTER 

The NASA SPoRT Center at the Marshall Space 
Flight Center (MSFC) seeks to accelerate the infusion of 
NASA Earth Science observations, data assimilation, 
and modeling research into weather forecast 
operations and decision-making at the regional and 
local level (Goodman et al. 2004). It directly supports 
the NASA strategic plan of using results of scientific 
discovery to directly benefit society. The program is 
executed in concert with other government, university, 
and private sector partners. The primary focus is on 
the regional scale and emphasizes forecast 
improvements on a time scale of 0-24 hours. The 
SPoRT Center has partnered with and facilitated the 
use of real-time NASA data and products to 17 
National Weather Service (NWS) Weather Forecast 
Offices (WFOs) primarily in the Southern Region, as 
well as several private weather entities. Numerous 
new techniques have been developed to transform 
satellite and lightning observations (Darden et al. 
2010) into useful parameters that better describe 
changing weather conditions, including proxy products 
that demonstrate utility for the upcoming GOES-R 
satellite era (Stano et al. 2010). 

The unique weather products have helped local 
WFOs improve forecasts of reduced visibility due to 
fog, low clouds, and smoke and haze from sources 
such as forest fires and agricultural burning, the onset 
of precipitation, the occurrence and location of severe 
weather events, and assess other local weather 
changes. Additionally, high-resolution satellite data 
provided by SPoRT has been used by the private sector 
to inform the marine weather community of changing 
ocean conditions and with tropical storm and 
hurricane monitoring. 

3. LAND INFORMATION SYSTEM DESCRIPTION 

The NASA LIS is a high performance land surface 
modeling and data assimilation system that integrates 
satellite-derived datasets, ground-based observations 
and model reanalyses to force a variety of LSMs 
(Kumar et al. 2006). By using scalable, high- 
performance computing and data management 
technologies, LIS can run LSMs offline globally with a 
grid spacing as fine as 1 km to characterize land 
surface states and fluxes. 

Case et al. (2008) presented improvements to 
simulated sea breezes and surface verification 
statistics over Florida by initializing the WRF model 
with land surface variables from an offline LIS spin-up 
run, conducted on the same WRF domain and 
resolution. In addition, Case et al. (2011) 
demonstrated the utility of using both the LIS land 
surface fields and high-resolution MODIS SSTs (Haines 
et al. 2007) to initialize the surface variables over the 
southeastern U.S., thereby providing a high-resolution 


lower boundary initial condition over the entire 
modeling domain that contributed to slight 
improvements in modeled summertime precipitation 
systems. 

To compare the SPoRT GVFs to the climatology 
GVFs, the LIS is configured to use the International 
Geosphere-Biosphere Programme (IGBP) land-use 
classification (Loveland et al. 2000) as applied to the 
MODIS instrument (Friedl et al. 2010). All static and 
dynamic land surface fields are masked based on the 
IGBP/MODIS land-use classes. The soil properties are 
represented by the State Soil Geographic (STATSGO; 
Miller and White 1998) database. Additional required 
parameters include quarterly climatologies of albedo 
(Briegleb et al. 1986), a 0.05° resolution maximum 
snow surface albedo derived from MODIS (Barlage et 
al. 2005), monthly climatologies of greenness fraction 
data derived from AVHRR in the control LIS simulations 
(Gutman and Ignatov 1998), and a deep soil 
temperature climatology (serving as a lower boundary 
condition for the soil layers) at 3 meters below ground, 
derived from 6 years of Global Data Analysis System 
(GDAS) 3-hourly averaged 2-m air temperatures using 
the method described in Chen and Dudhia (2001). 

4. EXPERIMENTAL NDVI/GVF DATASET 

4.1 MODIS NDVI Compositing algorithm 

The NASA SPoRT Center has been producing daily 
real-time, 0.01 degree (~l-km) resolution MODIS NDVI 
gridded composites over a Continental U.S. (CONUS) 
domain since 1 June 2010. The geographical extent of 
the CONUS grid ranges from 23°N to 52°N latitude and 
128°W to 65°W longitude. These composites are 
updated daily based on swath data from the MODIS 
sensor aboard the polar orbiting NASA Aqua and Terra 
satellites, with a product time lag of about one day. 

Figure 1 illustrates the data flow and processing 
that occurs during the course of creating a single day's 
NDVI composite. NDVI swath and cloud mask data are 
received from the University of Wisconsin, and each 
swath is individually mapped onto the CONUS 1-km 
grid using MclDAS utilities (left column of Figure 1). A 
simple time-weighting algorithm is applied to the 
remapped NDVI swath data (right column of Figure 1) 
that queries the previous 20 days of data to ensure a 
continuous grid of data populated at all pixels using 
the following formula at each (j, i) grid point: 
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NDVf(j,i) = \ 


H NDVI „ 
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;m > 1 


DaysLate n + l y 
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4.2 Calculating MODIS GVF for use in LIS 


The GVF from the MODIS NDVI data is calculated 
on the identical 0.01-deg grid following the procedures 
outlined in Zeng et al. (2000) and Miller et al. (2006). 
The GVF at each grid point is computed as a function of 
the IGBP/MODIS vegetation class by first determining 
the maximum NDVI at each grid point using the 
previous t NDVI composites, where t is ideally a 
running collection of the previous 365 days of NDVI 
composites. However, since a full year of NDVI 
composites has not yet been collected, we used all the 
composites from June to October to determine the 
NDVI max values at each grid point. In the future, the 
NDVI max values will be determined based on the 
previous year of daily NDVI composites. Once the 
NDVI max is found at each grid point, the NDVI max values 
are sorted as a function of land use class. All grid 
points with the same land use class are lumped 
together into a single distribution of NDVI max values, 
which are then sorted to find the 90 th percentile 
NDVI max for each land use class, and the 5 th percentile 
for the barren land use class. The GVF can then be 
computed as a function of vegetation (land use) class 
using the following formula: 


NDVI t - NDVI s 
NDVI Vi -NDVI S 


( 3 ), 


where NDVI , is the actual NDVI composite value at a 
grid point /, NDVI S is a global constant that corresponds 
to the 5 th percentile of the array of NDVI max values of 
the barren vegetation class, and NDVI V/i is the 90 th 
percentile of NDVI max values for the vegetation class at 
grid point /. 

Missing NDVI pixels are filled with the AVHRR 
monthly climatological GVF data, which should mainly 
impact winter composites when snow cover (or 
persistent cloud cover) prevents an adequate NDVI 
reading. Since evapotranspiration is much less 
substantial in the surface energy budget during the 
winter months, the impact of filling data with the 
AVHRR monthly climatology should be minimal. 

The new daily GVF dataset then replaces the 
monthly climatological GVF database based on five 
years of AVHRR observations currently available to the 
Noah LSM in both LIS and the public version of the 


WRF model. The much higher spatial resolution (1 km 
versus 0.15 degree) and daily updates based on real- 
time satellite observations have the capability to 
greatly improve the simulation of the surface energy 
budget in the Noah LSM within LIS and WRF. 

5. PRELIMINARY IMPACT RESULTS 

To measure the sensitivities and impacts, the LIS is 
configured to use the Noah LSM as run at NCEP and 
within the WRF model. For the initial tests, the LIS- 
Noah is run in an uncoupled, or offline mode with 
atmospheric analyses from the North American Land 
Data Assimilation System (NLDAS, Mitchell et al. 2004) 
and the Global Data Assimilation System (GDAS, 
Derber et al. 1991) providing the required input fields 
to drive the LSM integrations. The NLDAS is used 
where possible, while outside of the NLDAS domain, 
the GDAS analyses are invoked. The GDAS forcing 
fields of downward-directed longwave radiation, 
surface pressure, 2-m air temperature, and 2-m 
specific humidity are corrected topographically via 
lapse-rate and hypsometric adjustments using the 
elevation data differences between the LIS and native 
GDAS forcing grids (Cosgrove et al. 2003). 

A fine-scale model equilibrium state is not 
necessary for the current set of comparison 
simulations; therefore, a long-term spin-up run of LIS- 
Noah is not made at this time. Instead, the LIS is 
configured to run on a CONUS-scale 4-km grid identical 
to the real-time WRF grid run by the National Severe 
Storms Laboratory in support of Hazardous Weather 
Testbed and the Storm Prediction Center (Kain et al. 
2010). The CONUS 4-km offline LIS runs are made 
from 1 June to 31 October 2010, spanning the warm 
season months. In addition, a 1-km LIS run is made 
during the month of June centered on Montana, 
helping to depict the green-up throughout the month 
as well as the importance of higher-resolution GVF in 
areas of complex terrain. 

An example comparison of the June GVFs on the 
full-resolution 0.01° grid is given in Figure 2. Overall, 
the SPoRT MODIS GVFs tend to have higher greenness 
values across much of the High Plains, far southeastern 
U.S., Mexico, and the southwestern U.S. Areas that 
have lower GVFs in the SPoRT product include the 
Great Lakes and northeastern U.S., the Appalachians, 
and portions of the northern Rockies (Figure 2, bottom 
panel). Worth noting is the much greater amount of 
detail seen in the SPoRT GVF dataset, particularly in 
the inter-mountain West. The large-scale patterns are 
consistent, however, helping to affirm the soundness 
of the SPoRT GVFs. 

Zooming into the Montana high-resolution run 
reveals detailed differences between the AVHRR 
climatology and the SPoRT MODIS GVFs related to the 
ability to resolve complex terrain features. The AVHRR 
climatology appears quite smooth on this scale 
showing a minimum in GVF over the prairies of eastern 
Montana, northern Wyoming, western South Dakota, 
and central Idaho (Figure 3, top). A broad maximum 
exists along the Montana-ldaho border associated with 
the high terrain along the Continental Divide. 
Meanwhile, the SPoRT GVFs show much greater detail 



able to resolve small river basins with locally higher 
GVF, and relatively lower GVFs in the higher ridgetops 
in northwestern Wyoming and along the Montana- 
Idaho border (Figure 3, middle and bottom). The 
SPoRT GVFs generally depict much higher GVFs along 
the High Plains in the western Dakotas and eastern 
Montana (Figure 3, bottom). 

These variations in GVF directly impact the surface 
energy budget through the partitioning of sensible and 
latent heat fluxes. With the ability to better resolve 
the ridgetops, the LIS run with the SPoRT GVFs depict 
substantially lower latent heat fluxes up to 100 W m" 2 
over the high terrain where the GVFs are lower than 
the AVHRR climatology (Figure 4). Meanwhile, the 
latent heat fluxes are locally-regionally higher by about 
the same amount over parts of the High Plains where 
the SPoRT GVFs are higher (Figure 4, bottom). Worth 
noting in the difference field is the amount of local, 
detailed variations that are consistent with the terrain 
elevation features of this region (not shown). 

Finally, the higher GVFs over the High Plains of 
eastern Montana and the western Dakotas are 
consistent with fairly significant positive precipitation 
anomalies during May and June 2010. The top panel of 
Figure 5 shows that much of the High Plains 
experienced very substantial amounts of rainfall during 
May of 150% to as much as 600% of normal rainfall, 
especially northeastern Montana to western South 
Dakota. Meanwhile, western Montana and northern 
Idaho had generally normal to below normal rainfall 
during May where the SPoRT/MODIS GVFs are 
somewhat lower than the AVHRR climatology. The 
month of June 2010 saw much of the same, but with 
the entire region experiencing above-average 
precipitation (Figure 5, bottom). Precipitation is not 
necessarily the sole driver in producing anomalies in 
GVFs; however, on a regional scale there is certainly a 
strong qualitative correlation between the monthly 
precipitation anomalies and the deviations between 
the SPoRT/MODIS and AVHRR climatology GVFs. 

Additional preliminary results from the 4-km NSSL 
CONUS domain runs during the 2010 warm season will 
be shown in the companion presentation. 

6. SUMMARY AND FUTURE WORK 

Future efforts will involve incorporating the 
experimental SPoRT GVF data set into coupled model 
runs in which the LIS framework is called from within 
the Advanced Research WRF (ARW, Skamarock et al. 
2008) model. With LIS being tightly coupled to the 
ARW (Kumar et al. 2007), this framework enables a 
seamless incorporation of the experimental daily GVF 
data into the WRF model. Such experiments would 
demonstrate the impacts and utility of the high- 
resolution, daily updated vegetation data set by 
comparing the surface energy budget using the SPoRT 
GVF to that using the climatological means. In 
addition, validation of near-surface meteorological 
variables and quantitative precipitation will be 
conducted to quantify the level of improvement in the 
model forecasts. Also, the development of a global 1- 
km daily-updated MODIS GVF product would be 
beneficial to the operational weather community. 


Finally, SPoRT also seeks to collaborate with its 
partners and current users of the coupled LIS/ARW 
framework. The Air Force Weather Agency (AFWA) 
runs the LIS/ARW coupled framework operationally for 
the numerous regional model forecast domains across 
the globe. In fact, AFWA has been simultaneously 
experimenting with a similar such GVF data set in their 
coupled LIS/ARW model (Ruhge and Barlage 2011, this 
conference). It is the intent of SPoRT to develop and 
transition a high-quality dataset that can be used by its 
partner organizations to improve the short-term 
forecast process. 
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Figure 1. Diagram of the data and process flow of the SPoRT NDVI compositing algorithm. 
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Figure 2. Depiction of the greenness vegetation fraction (GVF) on the 
SPoRT CONUS 0.01° domain from 15 June for the current NCEP/AVHRR 
climatology (top), the SPoRT daily GVF (middle), and the difference 
(SPoRT-NCEP, bottom). 
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Figure 3. Depiction of the greenness vegetation fraction (GVF) on a 
0.01° resolution grid centered on Montana at 1800 UTC 27 June for 
the interpolated NCEP/AVHRR climatology (top), the SPoRT daily 
GVF (middle), and the difference (SPoRT - NCEP, bottom). 
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Figure 4. Depiction of the LIS-Noah latent heat flux (W m" 2 ) on a 0.01° 
resolution grid centered on Montana at 1800 UTC 27 June for the 
interpolated NCEP/AVHRR climatology (top left), the SPoRT daily GVF 
(top right), and the difference (SPoRT - NCEP, bottom panel). 
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Figure 5. Monthly percent of normal precipitation for May 2010 (top), and June 2010 (bottom), 
produced by the Advanced Hydrological Prediction Service and available at 
http://water.weather.gov/precip/ . 
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Importance of Vegetation in Models 


^^Evapotranspiration from Healthy Vegetation 

o Significant contribution of moisture transport into 
boundary layer 

o Important to represent vegetation density accurately in 
models 

Horizontal and Vertical Density of Vegetation 

o Greenness Vegetation Fraction (GVF), horizontal density 

o Leaf Area Index (LAI), vertical density 

o In Noah Land Surface Model (LSM), LAI is prescribed while 
GVF varies spatially by model grid cell 
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Land Information System (LIS) Overview 
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Current Operational GVF Dataset 


^^Five-Year Monthly Global Climatology 

o Derived from AVHRR NDVI data from 1985-1991 
o 0.144° resolution, valid at mid-point of each month 
o Currently in Noah LSM within NCEP/NAM and WRF models 
o Operational at NCEP since 1997 (Ek et al. 2003) 

Cannot account for variations in GVF due to: 

o Weather/climate anomalies (e.g. drought, excessive rain) 
o Land-use changes since the early 1990s (e.g. urban sprawl) 
o Wildfires and prescribed burn regions 

L J 
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SPoRT Real-Time Vegetation Products 


r 


Continental-U.S. NDVI grid at 0.01° resolution 



o Normalized Difference Vegetation Index (NDVI) from the 
Moderate Resolution Imaging Spectroradiometer (MODIS) 

o Daily swath data mapped to CONUS grid 

o NDVI composited using time-weighted formula and up to 
6 NDVI values in the previous 20 days: 
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SPoRT Real-Time Vegetation Products 


^Calculate GVF on 0.01° grid for use in LIS/Noah LSM^^ 

o Follows Zeng et al. (2000) and Miller et al. (2006) 


o Create distributions of NDVI max as a function of land use 
o Obtain 90 th percentile NDVI max for each land class (NDVI Vj ) 
o Determine 5 th percentile NDVI max for barren land use (NDVI S ) 
o Calculate GVF using the following formula: 


L 




NDVIyj - NDVI S 


o Implement into LIS software through a plug-in module 


J 
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Comparison of 15 June GVFs: AVHRR Climo 


NCEP/AVHRR June Climatology GVF (%) 
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Comparison of 15 June GVFs: SPoRT/MODIS 

SPoRT/MODIS GVF (%) on 15 June 2010 
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15 June GVF Difference (SPoRT - Climo) 


GVF Difference (SPoRT— NCEP) on 15 June 2010 



GVF Montana Zoom: 27 June 2010 
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LIS-Noah: Sensible HF at 18z 27 June 
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LIS-Noah: Latent HF at 18z 27 June 
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Animation of Daily GVFs: 1 June to 30 Oct 


SPoRT/UODIS GVF (%) valid 01 JUN 2010 
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Summary and Future Work 


r SPoRT Producing Real-time, Daily NDVI/GVFs 

o Run within the NASA LIS using the Noah LSM 
o Strong sensitivity of heat flux partitioning to GVFs 
o Can better resolve complex terrain than climo dataset 
o Adjusts realistically to anomalies in precipitation 

Future Efforts 

o Document improvements to fluxes, PBL evolution, & 
precipitation forecasts in coupled LIS-ARW runs 

o Examine detailed case studies through collaborations 

o Implement into real-time LIS-Noah/LIS-ARW runs 
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